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Abstract 



The results of a search for the production of second generation scalar leptoquarks are presented 
for final states consisting of either two muons and at least two jets or a muon plus missing transverse 
momentum and at least two jets. A total of 1.03 fb"^ integrated luminosity of proton-proton collision 
data produced by the Large Hadron Collider 3X ^/s = 1 TeV and recorded by the ATLAS detector is 
used for the search. The event yields in the signal regions are found to be consistent with the Standard 
Model background expectations. The production of second generation leptoquarks is excluded for a 
leptoquark mass m-LQ < 594 (685) GeV at 95% confidence level, for a branching ratio of 0.5 (1 .0) for 
leptoquark decay to a muon and a quark. 
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Abstract The results of a search for the production 
of second generation scalar leptoquarks are presented 
for final states consisting of either two muons and at 
least two jets or a muon plus missing transverse mo- 
mentum and at least two jets. A total of 1.03 fb~^ inte- 
grated luminosity of proton-proton collision data pro- 
duced by the Large Hadron Collider at y/s = 7 TeV 
and recorded by the ATLAS detector is used for the 
search. The event yields in the signal regions are found 
to be consistent with the Standard Model background 
expectations. The production of second generation lep- 
toquarks is excluded for a leptoquark mass tolq < 594 
(685) GeV at 95% confidence level, for a branching ra- 
tio of 0.5 (1.0) for leptoquark decay to a muon and a 
quark. 



1 Introduction 

The remarkable similarities between quarks and leptons 
in the Standard Model (SM) lead to the supposition 
that there could be a fundamental relationship between 
them at a sufficiently high energy scale, manifested by 
the existence of leptoquarks (LQ) [T]. LQs are hypothet- 
ical particles which carry both baryon and lepton num- 
ber and have fractional electrical charge. The present 
search is performed within the minimal Buchmiiller- 
Riickl-Wyler model 2 , where LQs are restricted to 
couple to quarks and leptons of one generation. In this 
model, LQs are required to have pure chiral couplings 
to SM fermions in order to avoid inducing four-fermion 
interactions that would cause flavour-changing neutral 
currents and lepton family-number violations. At the 
Large Hadron Collider (LHC), scalar LQs can be pro- 



duced either in pairs or singly. Single LQ production 
involves the unknown X^g-i-q coupling, while pair pro- 
duction of scalar LQs occurs mostly via gluon-gluon fu- 
sion, dominant for tolq ^ 1 TeV, and gg-annihilation, 
dominant at larger masses. Both pair-production modes 
involve only the strong coupling constant, and therefore 
all model dependence is contained in the assumed LQ 
mass JTiLQ and the branching ratio /3 for LQ decay to a 
charged lepton and a quark [j. LQs can also decay to a 
neutrino and a quark; in this case, the branching ratio 
is 1 — /3. Pair production of scalar LQs at the LHC has 
been calculated at next-to-leading order (NLO) \£^. 

The results presented in this paper are an update 
of the previous ATLAS search for second generation 
LQs [5] and extend the bounds arising from previous di- 
rect searches performed by CMS ^, ATLAS 0, DO [7] 
and OPAL [5]. A total integrated luminosity of 1.03 fb~^ 
of proton-proton collision data at a centre of mass en- 
ergy y/s = 7 TeV, collected with the ATLAS detector 
from March through July 2011, is used for the search. 
The final states arising from leptoquark pairs decay- 
ing into two muons and two quarks (fJ-fJ-jj), or into a 
muon, a neutrino and two quarks {fJ-vjj), are consid- 
ered. These result in experimental signatures of either 
two high transverse momentum (pt) muons and two 
high pt jets, or one high px muon, missing transverse 
momentum, and two high px jets. 

Analyses for both dimuon and single muon final 
states start with the selection of event samples with 

^ The \i^Q_i_q coupling determines the LQ lifetime and 
width [3]. For LQ masses considered here, 200 GeV< 
mLQ <700 GeV, couplings greater than e x 10~®, with e = 
•/Ana the electron charge, and a{Mz) = 1/128, correspond 
to decay lengths less than roughly 1 mm. In addition, to be 
insensitive to the coupling, the width cannot be larger than 
the experimental resolution of a few GeV. This sets the ap- 
proximate sensitivity to the unknown coupling strength. 
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large signal acceptance. Since background cross sec- 
tions are several orders of magnitude larger than the 
signal cross sections, these samples are dominated by 
the major backgrounds: Z+jets and tt in the /^/ijj case, 
and Ty+jets and tt for the fii'jj case. Further selec- 
tion requirements are then applied to these samples to 
define control regions used to determine the normal- 
ization of the aforementioned backgrounds. The deter- 
mination of the multi-jet background is performed in 
a fully data-driven approach, and the smaller diboson 
and single top-quark backgrounds are estimated using 
Monte Carlo (MC) simulations. 

After all background contributions are determined, 
variables selected to enhance the discrimination between 
signal and background are combined into a log like- 
lihood ratio, which is used to search for an excess of 
events over the SM background prediction. The searches 
are performed independently for each final state. The 
results are then combined and interpreted as lower bounds 
on the LQ mass for different /? hypotheses. 



tum measurements. The niuon spectrometer is based 
on three large superconducting toroids with coils ar- 
ranged in an eight-fold symmetry around the calorime- 
ters, covering a range of |?7| < 2.7. Over most of the rj 
range, precision measurements of the track coordinates 
in the principal bending direction of the magnetic field 
are provided by Monitored Drift Tubes (MDTs). At 
large pseudorapidities (2.0 < \ri\ < 2.7), Cathode Strip 
Chambers (CSCs) with higher granularity are used in 
the innermost station. 

A three-level trigger system selects events to be recorded 
for offline analysis. The muon trigger detectors consist 
of Resistive Plate Chambers (RPCs) in the barrel (|?7| < 
1.05) and Thin Gap Chambers (TGCs) in the end-cap 
regions (1.05 < \r]\ < 2.4), with a small overlap in the 
JT^I = 1.05 region. The data considered in this analysis 
are selected from events containing at least one muon 
with the transverse momentum determined by the trig- 
ger system satisfying pT > 18 GeV. 



2 The ATLAS detector 

The ATLAS detector [S] is a multi-purpose detector 
with a forward-backward symmetric cylindrical geome- 
try and nearly 47r coverage in solid angle q 

The three major sub-components of ATLAS are the 
tracking detectors, the calorimeters and the muon spec- 
trometer. Charged particle tracks and vertices are re- 
constructed with silicon-based tracking detectors that 
cover |ry| < 2.5 and a transition radiation tracker ex- 
tending to \ri\ < 2.0. The inner tracking system is im- 
mersed in a homogeneous 2 T axial magnetic field pro- 
vided by a solenoid. Electron, photon, and jet energies 
are measured in the calorimeters. The calorimeter sys- 
tem is segmented into a central barrel and two endcaps, 
collectively covering the pseudorapidity range of |?7| < 
4.9. A liquid-argon (LAr) electromagnetic calorimeter 
covers the range JTyj < 3.2 and an iron-scintillator tile 
hadronic calorimeter covers the range |?7| < 1.7. Endcap 
and forward LAr calorimeters provide both electromag- 
netic and hadronic measurements and cover the region 
1.5 < |?7| <4.9. 

Surrounding the calorimeters, a muon spectrome- 
ter [5] with air-core toroids, a system of precision track- 
ing chambers, and detectors with triggering capabili- 
ties provides muon identification and precise momen- 

^ ATLAS uses a right-handed coordinate system with its 
origin at the nominal interaction point and the z-axis along 
the beam pipe. Cylindrical coordinates {r,(f>) are used in the 
transverse plane, with ^ the azimuthal angle around the beam 
pipe. The pseudorapidity r? is defined in terms of the polar 
angle 6 hy r] = -In tan(S/2). 



3 Simulated samples 

Simulated event samples are used to determine all sig- 
nal and some of the background yields. Signal samples 
for LQ masses between 200 GeV and 1000 GeV are sim- 
ulated with PYTHIA 6.4.25 [TO]. NLO cross sections as 
determined in Ref. H], using CTEQ6.6 [TT] parton dis- 
tribution functions (PDFs), are used to normalize the 
samples at each mass point. 

Samples of W and Zf-f* production in association 
with n partons (where n can be 0, 1, 2, 3, 4 and 5 
or more) are simulated with the ALPGEN [T^] gen- 
erator interfaced to HERWIG [H] and JIMMY [H] to 
model parton showers and multiple parton interactions, 
respectively. The MLM 12 parton-shower matching 
scheme is used to form inclusive W^/Z-|-jets samples. 
MCQNLO [TS] is used to estimate the production of 
single top quarks and top quark pairs. A top quark mass 
of 172.5 GeV is used in the simulation. Diboson events 
are generated using HERWIG, and the cross sections 
are scaled to NLO calculations |15[|TB] . 

All simulated events are passed through a full detec- 
tor simulation based on GEANT4 [T7] and then recon- 
structed with the same software chain as the data |18] . 
During the data-taking period considered in this search, 
the mean number of primary proton-proton interactions 
per bunch crossing was approximately six. The effect of 
this pile-up is taken into account in the analysis by over- 
laying simulated minimum bias events onto the simu- 
lated hard-scattering events. The MC samples are then 
reweighted such that the average number of pile-up in- 
teractions matches that seen in the data. 
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4 Object and Event selection 

Collision events are identified by requiring at least one 
reconstructed primary vertex candidate with at least 
three associated tracks with px. track > 0.4 GeV. If two 
or more such vertices are found, the one with the largest 
sum oip^ track ^^ taken to be the primary vertex. Muons 
are reconstructed by matching tracks in the inner de- 
tector to track segments in the muon spectrometers, as 
described in Ref. [19 . In addition to the track qual- 
ity requirements imposed for identification, the muon 
tracks must also satisfy jdoj < 0.1 mm and |zo| < 5 mm, 
where do and zq are the transverse and longitudinal im- 
pact parameters measured with respect to the primary 
vertex. All selected muons must have pr > 30 GeV and 
are restricted to be within \ri\ < 2.4. Muon candidates 
must pass the isolation requirement p™"°^°/pT < 0.2, 
where p^°"°^" is the sum of the px of the tracks within 
AR = a/(Z\(/))2 + {Af]'^) < 0.2 of the muon track, ex- 
cluding the muon px contribution. Selected events must 
have at least one muon identified by the trigger system 
within a cone AR < 0.1 centered on a selected muon. 

Jets are reconstructed from calorimeter energy clus- 
ters using the anti-fci algorithm [2T\ with a radius pa- 
rameter R = 0.4. Corrections are applied in order to ac- 
count for the effects of the non-compensating calorime- 
ter, upstream material and other effects, by using px 
and ry-dependent correction factors derived from simu- 
lation and validated with test-beam [22] and collision 
data studies 123,. After applying quality requirements 
based on shower shape and signal timing with respect 
to the beam crossing, the selected jets must satisfy 
Pt > 30 GeV, 1 77 1 < 2.8 and must be separated from the 
selected candidate muons by AR ^0.4. The presence 
of neutrinos is inferred from the missing transverse mo- 
mentum i?™'^^, defined as the magnitude of the negative 
vector sum of the transverse momenta of reconstructed 
electrons, muons and jets, as well as calorimeter energy 
deposits not associated to reconstructed objects. 

Corrections to the muon trigger and reconstruction 
efficiencies and to the momentum resolution are applied 
to the simulated events so that their kinematic distribu- 
tions match those observed in data, with an impact on 
the predicted number of events of less than 2%. These 
corrections are derived from samples of Z — > fifj, and 
W ^ fiv decays JS.;, taking into account the effects of 
multiple scattering and the intrinsic resolution of the 
muon spectrometer [20] . In order to validate the cor- 
rections at high px, the alignment of the muon spec- 
trometer, which dominates the momentum resolution 
for px larger than approximately 200 GeV, is derived 
from a sample of straight track data taken in special 



runs with the toroids turned off, resulting in agreement 
within the considered systematic uncertainties. 

Events selected for this search are required to con- 
tain either exactly two muons and at least two jets for 
the /i/ijj final state, or exactly one muon, at least two 
jets and E^^^^ > 30 GeV for the ^vjj final state. In 
the /i/ijj channel, only events with m^^ > 40 GeV are 
considered. In the ^vjj channel, events are required to 
have TOx = ^j2pi^E^''"'(l - cos(Z\</>)) > 40 GeV, where 
Acf) is the angle between the muon and the i?™'^** di- 
rection in the plane perpendicular to the beam. Events 
with identified electrons as defined in Ref. [24], with 
px > 30 GeV, and I77I < 2.47 are rejected. After all 
the selection criteria are applied the acceptance times 
efficiency ranges from about 60% (55%) for a LQ signal 
of TOLQ = 300 GeV to 65% (60%) for a LQ signal of 
wlq = 600 GeV for the ^^.jj (fJ-i^jj) channel. 



5 Background determination 

Major backgrounds in this search arise from y-|-jets 
(y = W, Z) and tt processes. The kinematic distribu- 
tions of these are determined using MC samples, and 
their absolute normalization is evaluated from data us- 
ing control regions, which are subsets of the selected 
sample, designed to enhance either the T^-l-jets or the 
top quark contribution. The multi-jet background is ob- 
tained directly from data and prior to the estimation of 
the normalization for the two main backgrounds, while 
the determination of the remaining backgrounds (dibo- 
son and single top quark production) relies entirely on 
MC simulations. 

Two control regions are used in the /Jfijj channel. 
(I) Z-fjets: formed by events within a narrow dimuon 
invariant mass m^^ window around the Z boson mass, 
defined by 81 < ni^j^j < 101 GeV, and at least two jets, 
and (II) tt: one of the muons is replaced by an electron 
resulting in events with a muon and an electron, and at 
least two jets. 

Three control regions are used in the iii'jj chan- 
nel. (I) W+ 2 jets: events in the vicinity of the W bo- 
son Jacobian peak, selected by requiring 40 < tox < 
120 GeV, exactly two jets and St < 225 GeV, where 
S'x is the scalar summed transverse energy ^x, defined 
as St = Pt + -E^'"" + P^x" + Px*^ (H) ^+ 3 jets: 
events passing the 40 < mx < 120 GeV requirement, 
with at least three jets and St < 225 GeV, and (III) 
tt: events with at least four jets, with pl^ > 50 GeV 
and pl^ > 40 GeV. In all of the control regions the 
expected signal yields are negligible. 

The normalizations of the y-|-jets and ii backgrounds 
are obtained by comparing data and MC yields in the 



The ATLAS Collaboration 



o ^ 

05 

10 

1 

10-^ 



ATI AC , ^r^ — - Data (\/s=7 TeV) 

ATLAS LQLQ^^,^, nv+Jets 

1 .03 fb"^ ■ Top 

I I Multi-jet 
■ Diboson 
— LQ (m=600 GeV) 



JLdt 




100 



200 



300 



400 500 600 
m(^,M ) [GeV] 



ATLAS LQLQ^nnJj 



^i- Data (\/s=7 TeV) 
nv+Jets 



Ldt = 1 .03 fb □ Top 




Multi-jet 
Diboson 
LQ (m=600 GeV) 



100 200 300 400 500 600 700 800 900 
Average LQ Mass [GeV] 



(a) 



(b) 



0)1' 



ATI AC . ^r^ -- Data (\/s=7 TeV) 

ATLAS LQLQ^mi □v.Jets 
Ldt = 1.03 fb"^ □ Top 
S^ -' Multi-jet 

"% ■ Diboson 

'H^ — LQ (m=600 GeV) n 




200 400 600 800 1000 1200 1400 1600 1800 

St [GeV] 



(c) 



Fig. 1 Distributions of the input LLR variables for the fif^jj channel for data and the SM backgrounds. I(a)| Invariant mass of 
the two muons in the event, [(b)] Average LQ mass resulting from the best muon-jet combinations in each event, and |(c)| 5T. 
The stacked distributions show? the various background contributions, and data are indicated by the points with error bars. 
The 600 GeV LQ signal is also shown for /3 = 1.0. In all figures, the last bin contains the sum of all entries equal to and above 
the bin lower boundary. 



control samples defined above. In the fJ-fxjj channel, 
each correction factor is obtained independently for each 
background, on account of the high purity of the two 
different control regions. In the fJ-i^jj channel, there is 
significant cross-region contamination and therefore the 
number of y+jets and ti events is determined by simul- 
taneously minimizing the x^ formed by the differences 
between the observed and predicted SM yields in the 
three control regions. The resulting scale factors are of 
the order of 10% in the low St region. 

The multi-jet background in the selected sample and 
in each control sample is obtained from a fit to the ttt,^^ 
and E^^^'^ distribution in the fj.fijj and fJ-i'jj channels, 
respectively. In these fits, the relative fraction of the 
multi-jet background is a free parameter, and the sum 
of the total predicted events is constrained to be equal 



to the total observed number of events. The y-|-jets and 
tt normalizations are not fixed. Multi-jet background 
arises predominantly from muons from secondary de- 
cays. Therefore, templates for the multi-jet background 
distributions are constructed from multi-jet enhanced 
samples of data events in which the muons fail the re- 
quirement on the transverse impact parameter or the 
isolation selection requirements described in Section SI 
In the fifj.jj channel, the W^-|-jets contribution is esti- 
mated together with the multi-jet background. During 
this procedure, the V^-l-jets and tt normalizations are 
fitted as well, providing an independent estimate. The 
resulting values agree with those obtained from the con- 
trol regions, which are the ones used in the analysis. 

After analyzing 1.03 fb~^ of data and applying the 
analysis requirements described in Sectional good agree- 
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Fig. 2 Distributions of the input LLR variables for the fj,ujj 
of the muon and the E!^^^^ in the event , [(b)] St , [(c)] LQ mass, 
various background contributions, and data are indicated by 
LQ signal is also showfu for /3 = 0.5. In all figures, the last bin 
boundary. 



channel for data and the SM backgrounds, [(a)] Transverse mass 
and 1(d) I LQ transverse mass. The stacked distributions show the 
the points with error bars. The expected signal for a 600 GeV 
contains the sum of all entries equal to and above the bin lower 



ment is observed between the data and the SM expecta- 
tion. The observed and expected yields in the selected 
sample are 9254 and 9300 ± 1700 for the /i/J.j j channel, 
and 97113 and 97000 ± 19000 for the fii^jj channel. For 
a LQ mass of 600 GeV, 8.2 ±0.4 and 3.9±0.2 events are 
expected for the fJ.fJ.jj and the fJi'jj final states, respec- 
tively. The aforementioned uncertainites fully account 
for (the dominant) systematic and statistical uncertain- 
ties. 



6 Likelihood analysis 

Several kinematic variables, selected to provide the best 
discrimination between LQ events and SM backgrounds, 
are combined in a log likelihood ratio in order to search 
for a LQ signal. In the fifjjj channel, m^^, St 
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jetl 

Ft 



Ml 

Pt 



pl^ and the average reconstructed lep- 



toquark mass rhhQ are used. In the fJiyjj channel, St, 
rriT, the transverse leptoquark mass TTirp^ and the lep- 
toquark mass m-LQ are used. The distributions of these 
input variables are shown in Fig. [T] and Fig. [2] for the 
fjfjjj and the fJiyjj final states, respectively. 

In the fjfjjj channel, an average LQ mass jtilq is de- 
fined for each event by reconstructing all possible com- 
binations of lepton-jet pairs, using the two highest pt 
jets in each event. Of the four possible combinations in 
each event, the pairing which provides the smallest dif- 
ference between the LQ masses is chosen, and their av- 
erage is used in the likelihood analysis. In the fJiyjj final 
state, because the longitudinal component of the neu- 
trino momentum is unknown, only one mass from the 
muon and a jet can be reconstructed, and the E!f^^^ and 
the remaining jet are used to calculate the transverse 
mass of the other LQ. The two masses which provide 
the smallest absolute difference are used in the likeli- 
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Table 1 The predicted and observed yields and the expected 
yields for a LQ signal of "Ilq = 600 GeV after requiring 
LLR > 2 for the fifijj channel and LLR > 7 for the f^ujj 
channel. The ^i^ijj {lJ.i^jj) channel signal yields are computed 
assuming 13 = 1.0 (0.5). Statistical and systematic uncertain- 
ties as described in Section[7]are shown. These are calculated 
assuming a 100% correlation for the same source between 
the different backgrounds. These systematic uncertainties are 
computed as the sum of the absolute values of the systematic 
variation in each bin and are shown to indicate the scale. This 
is an approximation to the standard ensemble method used 
in the limit setting code. 

Source fJ.fJ.jj Channel fJi^jj Channel 



y+jets 
Top 

Diboson 
Multi-jet 



14.2 ±6.4 

3.0 ±2.2 

0.8 ±0.6 

<0.1 



12.9 ±9.9 

1.9 ±1.2 

0.3 ±0.1 

<0.1 



Total 
Data 



18 ±8 
16 



15 ±11 
14 



LQ 



8.2 ±0.4 



3.2 ±0.2 



hood analysis. With this algorithm, the probability of 
picking the correct pairing is of around 90% for both 
channels. 

For each event, likelihoods are constructed for the 
background (Lb) and the various signal LQ hypothesis 
(Ls) as follows: Lb = Hb^ixij), Ls = 11 Si{xij), where 
hi, Si are the probabilities of the i-th input variable 
from the normalized summed background and signal 
distributions, respectively, and Xij is the value of that 
variable for the j'-th event in a sample. The log likeli- 
hood ratio for each tested signal, LLR —log{Ls/LB), 
is used as the final variable to search for the LQ signal. 



7 Systematic uncertainties 

Systematic uncertainties originating from several sources 
are considered. These include uncertainties in lepton 
momentum, jet energy and E^^^^ scales and resolutions 
and their dependence on the number of pile-up events, 
the background estimations, and the LQ production 
cross section. For each source of uncertainty consid- 
ered, the analysis is repeated with the relevant variable 
varied within its uncertainty, and a new LLR is built 
for the systematically varied sample, enabling the un- 
certainty in both the predicted yield and the kinematic 
distributions to be propagated to the final result. In this 
section, systematic uncertainties are described for each 
source of systematics, calculated assuming each source 
to be 100% correlated among the different backgrounds. 
Uncertainties are given for the region of LLR > 2 and 
LLR > 7 for the /Lt/zjj and the iJivjj channels, respec- 
tively, although the full LLR distribution is used to 
search for the LQ signal. 
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Fig. 3 I (a) I LLR distributions for the tJ.tJ.jj and |(b)| for the 
fjiyjj final states for a LQ mass of 600 GeV. The data are 
indicated with the points and the filled histograms show the 
SM background. The multi-jet background is estimated from 
data, while the other background contributions are obtained 
from simulated samples as described in the text. The LQ 
signal corresponding to a LQ mass of 600 GeV is indicated 
by a solid line, and is normalized assuming /3 = 1.0(0.5) in 
the fJfJJJ ifJ'^Jj) channel. The lowest bin corresponds to back- 
ground events in regions of the phase space for which no signal 
events are expected. 



The jet energy scale ( JES) and resolution ( JER) are 
varied up and down by Ict [23] for all simulated events. 
Their impact is estimated independently, and the cor- 
responding variations are propagated to the E!^^^'^ in 
the case of the iJ.iyjj channel. The resulting effect of the 
JES (JER) uncertainty is 9% (8%) and 15% (7%) for 
the backgrounds in the ii^ijj and the iivjj channels, 
respectively. For a LQ signal of tolq = 600 GeV, both 
are 1% for the hijLJj channel, and 2.4% and 1% for the 
^lyjj channel. 

The systematic uncertainties from the muon reso- 
lution and momentum scale are derived by comparing 
the m^^ distribution in Z =>■ /i/j, control samples to 
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Z —i' fifj, MC samples and are approximately 1% pU] . 
These result in uncertainties of 12% and 3% for the 
total background prediction in the fj.fj.jj and the fjvjj 
channels, respectively, and in uncertainties of 1.4% for 
a LQ signal of tolq = 600 GeV for the fifijj and the 
fivjj channels. 

Systematic uncertainties due to assumptions in the 
modelling of the y+jets background are estimated by 
using SHERPA [IS] samples instead of the ALPGEN 
samples described in Section |31 The resulting uncer- 
tainty is 30% for the fjfjjj channel and 60% for the 
fj.iyjj channel. Similarly, systematic uncertainties aris- 
ing from the modelling of the tt process are obtained by 
using different parameter values to simulate alternative 
samples to the one described in Section [3] These in- 
clude samples in which the top quark mass is varied up 
and down by 2.5 GeV, generated with MC@NLO, sam- 
ples where the initial and final-state radiation (ISR and 
FSR) contributions are varied accordingly to their un- 
certainties, generated with ACER MC [2S] , and samples 
generated with POWHEG [17] interfaced to PYTHIA 
and JIMMY. These impact the total background yields 
by 12% (7%) for the fifijj {j^vjj) final state. For both 
V-hjets and tt backgrounds, a 10% uncertainty on the 
scale factors is considered, covering the variation of the 
scale factors in the low and high px regions. 

Systematic uncertainties in the multi-jet background 
in the fjfijj channel are determined by comparing re- 
sults derived from fits to kinematic variables other than 
the nominal ones. These include the leading muon px, 
the leading jet px, the E^^^^ and the scalar sum of 
the transverse momenta of the two muons in the event. 
In the fJivjj channel, an alternative loose-tight matrix 
method [5^ with two different multi-jet enhanced sam- 
ples obtained by inverting the isolation and the |do| 
requirements is used. Since the relevant phase space of 
the multi-jets in the two channels is very different, the 
different control regions have very different statistics 
which leads to a large difference in precision to which 
this background can be estimated. The resulting uncer- 
tainties are 90% in the fifiji channel and 33% in the 
fivjj channel. 

A luminosity uncertainty of 3.7% [H] is assigned 
to the LQ signal yields and to the yields of background 
processes determined from simulation: diboson and sin- 
gle top quark production. Further systematic uncertain- 
ties considered arise from the finite number of events in 
the simulated samples, amounting to 4%-25% depend- 
ing on the LQ mass being considered. 

For the signal samples, additional systematic uncer- 
tainties originate from ISR and FSR effects, resulting in 
an uncertainty of 2% for both channels. The choice of 
the renormalization and factorization scales, which are 



varied from ttt-lq to 2tolq and ?7^LQ/2, and the choice 
of the PDF, determined with the CTEQ eigenvectors 
errors and by using the MRST2007LO* PDF set [30] . 
result in an uncertainty in the signal acceptance of 1%- 
6% for LQ masses between 300 GeV and 700 GeV. 



8 Results 

Figure [3] shows the LLR for the data, the predicted 
backgrounds and a LQ signal of 600 GeV for the fifij] 
and the fJvjj channels. To ensure sufficient background 
statistics, bins with a total background yield less than 
twice the statistical uncertainty in that bin are merged 
into a single bin. There is no significant excess in data 
observed at large LLR values where such a signal would 
appear, and the data are found to be consistent with 
the SM background expectations (see Table [Ij . Upper 
limits are derived at 95% confidence level (CL) for the 
scalar leptoquark production cross section using a mod- 
ified frequentist CLs approach '311132]. The test statis- 
tic is defined as — 21n(Q) = — 21n(Ls+{,/L{,), where the 
likelihoods Ls+fc and L^ follow a Poisson distribution 
and are calculated based on the corresponding LLR 
distributions. Systematic uncertainties as described in 
Section [7] are treated as nuisance parameters with a 
Gaussian probability density function. 

The 95% CL upper bounds on the cross section for 
leptoquark pair production as a function of mass are 
shown in Fig. 2] for the fifijj and the fivjj channels at 
/3 = 1.0 and (5 = 0.5, respectively. The expected and 
observed limits for the combined channels are shown in 
the P vs. rTiLQ plane in Fig. [5] 



9 Conclusions 

The results of a search for the pair production of sec- 
ond generation scalar leptoquarks using 1.03 fb^^ of 
proton-proton collision data produced by the LHC at 
-y/s = 7 TeV and recorded by the ATLAS detector are 
presented. The data are in good agreement with the 
expected SM background, and no evidence of LQ pro- 
duction is observed. Lower limits on leptoquark masses 
of TOLQ > 685 GeV and mLQ > 594 GeV for ^ = 1.0 
and /3 — 0.5 are obtained at 95% CL, whereas the ex- 
pected limits are tolq > 671 GeV and ttilq > 605 GeV, 
respectively. These are the most stringent limits to date 
arising from direct searches for second generation scalar 
leptoquarks. 
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Fig. 4 |(a)| 95% CL upper limit on the pair production cross 
section of second generation leptoquarks for the /^/ijj channel 
at fS =1.0 and | (b) | for the fj.iyjj channel at /3 =0.5. The solid 
lines indicate the individual observed limits, while the ex- 
pected limits are indicated by the dashed lines. The theoret- 
ical prediction is indicated by the hatched band and includes 
the systematic uncertainties due to the choices of the PDF 
and the renormalization and factorization scales. The dark 
(green) and light (yellow) solid band contains 68% (95%), re- 
spectively, of possible outcomes from pseudo-experiments in 
which the yield is Poisson-fluctuated around the background- 
only expectation. 
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